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We reinvestigate the large degeneracy solution of the multichannel Kondo problem, and show 
how in the universal regime the complicated integral equations simplifying the problem can be 
mapped onto a first order differential equation. This leads to an explicit expression for the full 
zero-temperature scaling functions at - and away from - the intermediate Non Fermi Liquid fixed 
point, providing complete analytic information on the universal low - and intermediate - energy 
properties of the model. These results also apply to the widely-used Non Crossing Approximation 
of the Anderson model, taken in the Kondo regime. An extension of this formalism for studying 
finite temperature effects is also proposed and offers a simple approach to solve models of strongly 
correlated electrons with relevance to the physics of heavy fermion compounds. 



Models of quantum impurities constitute a central 
piece in our understanding of strongly correlated systems, 
as demonstrated by the importance of these ideas for 
the physical description of magnetic impurities in met- 
als or superconductors 0jl3, the correlation-driven 
metal to insulator transition the behavior of heavy 
fermion compounds 0, Ifl, and the coherent transport 
through nanostructures pj ■ Despite of their relative sim- 
plicity, one usually has recourse to different techniques 
in order to fully grasp the complete physics displayed by 
these localized hamiltonians. Unfortunately, analytical 
approximation schemes devised to deal with the Kondo 
problem, like the slave boson method @, lead in gen- 
eral to too much simplifications. On the contrary, purely 
numerical methods Q or exact techniques [To| provide 
quantitative results, but often lack the flexibility to solve 
more interesting cases, like lattice impurity models, or to 
compute dynamical observables. Conformal field theory 
(CFT) can also provide interesting insights on this prob- 
lem 1 1 1| . which is however limited to the elucidation of 
the critical behavior right at the possible fixed points. 

Another way of tackling impurity hamiltonians lies in 
setting up controlled limits that lead to a simplified so- 
lution encompassing all physical aspects of the problem. 
We will focus here on a large number of channels limit 
of quantum impurity models that was pro pos ed p revi- 

ously mm, 

from which the well-studied [lj, [15| self- 
consistent pertubation theory (Non Crossing Approxima- 
tion or NCA) appears to be a limiting case. However sim- 
plified, this alternative approach is known to lead to chal- 
lenging non-linear integral equations, from which limited 
analytical information can be gathered |12lll6j . A numer- 
ical solution of these equations is often possible, but is 
increasingly difficult to obtain at low temperature and for 
real frequency quantities. Thus little work has been un- 
dertaken using this method for the many interesting ex- 
tensions of the Kondo model, like multi-impurity hamil- 
tonians or the Kondo lattice. This is unfortunate, as 
this technique offers a concrete opportunity for building a 
consistent theory for the competition of Kondo screening 
and magnetism mi which lies at the heart of the quan- 



tum criticality displayed in heavy fermions compounds. 

We will show in this paper that much progress can be 
made towards both analytic and numerical solutions of 
the multichannel Kondo model. Indeed we are able for 
the first time to solve this problem at zero temperature in 
the universal regime (when the Fermi energy is large com- 
pared to frequency). The key step of this solution lies in 
the mapping of the saddle-point integral equations onto a 
much simpler non-linear differential equation, that can be 
solved explicitely. Our calculation follows the well-known 
analysis by Mueller-Hartmann of the similar NCA equa- 
tions for the infinite U Anderson impurity model |16| . 
However the Kondo regime of these equations seems not 
to have been considered previously in the literature, and 
constitutes the originality of the present work with re- 
spect to mi- We are thus able to obtain analytic ex- 
pressions for the frequency-dependent scaling functions, 
for the Kondo temperature and for the precise location 
of the intermediate Non Fermi Liquid fixed point. This 
offers also new insights on quantum criticality: we intro- 
duce below a concept of complex Kondo temperature and 
describe an interesting crossover phase diagram around 
the intermediate coupling fixed point. We will finally 
show how these ideas can be formulated at finite temper- 
ature, by deriving non-linear finite difference equations 
for the Green's functions, that may be useful to solve 
numerically more challenging extensions of the present 
problem. 

Considering the Kondo model from now on, we will be 
interested here in a multichannel version that consists of a 
SU(iV) spin antiferromagnetically coupled to K channels 
of conduction electrons, as expressed by the hamiltonian: 

H=Y^ £kcl aa c kaa + S a , a c\ aa c k , ala (1) 

k,a,a kk' aa' a 

using a fermionic representation S a r a = f\,fa — qo5 aa i 
of the localized spin, where a = 1 . . . N . The constraint 
Ta fa fa = QoN should also be enforced. Here the pa- 
rameter qoN may be understood as the size of the quan- 
tum spin, and the channel index a is taken to run from 
1 to K. 
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In a previous work [12|], Parcollet and Georges have 
shown how this model could be simplified in the limit 
where both N and K are large, keeping 7 = K/N fixed 
(the fact that the spin size qoN is taken to be large is also 
an important point in obtaining their solution). Intro- 
ducing the spinon propagator G/(r) = — (T T / CT (r)/ ( J(0)) 
and the Green's function G b (t) = +(T t B q (t)BI(0)) of 
the channel bosonic field B a (r) conjugate to the operator 
Sfco- c ko-af<7> t ne saddle-point equations read (in imagi- 
nary time): 

E/(r) = 7 Go(r)G B (r) (2) 
£ B (r) = -G (-T)G f (r) = G (T)G f (r) (3) 

for a particle-hole symmetric bath Green's function 
Go(iuj) = Y^k — £fc)- The self-energies in the previ- 

ous set of equations are given by the Dyson equations: 

Gj (iw n ) = iu n + A - S/(iw n ) (4) 

Gg'tifn) = T, B (iv n ) (5) 

where oj n (y n ) denotes a fermionic (bosonic) Matsubara 
frequency. The Lagrange multiplier A is used to enforce 
the constraint related to the spin size, which reads in the 
large N limit: 



G f (T = Q-) = q 



(6) 



The usual NCA equations 1_4( are obtained when the 
constraint flf a = 1 is performed exactly, and this 
amounts to take the qo — ► limit in the saddle-point 
equations. It is useful to analytically continue equa- 
tions J5J- @, and one readily obtains in terms of retarded 
quantities at zero temperature [l7| : 



"+00 j 

-/ - G '^e)[G B {u-e)-J K ) 
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7T 

+00 1 

de 



G (e)G B (c-e) 



^G^G^-e) 



+ °° de 

-G (e)G';( W - e) 



(7) 



(8) 



where double primed quantities denote imaginary parts. 
The subtracted Jk in comes from a careful inspec- 
tion of the short time behavior of G b (t). As we are 
interested in the universal regime where the Fermi en- 
ergy of the conduction electrons is the biggest energy 
scale, we follow and consider a flat density of states 
of half-width A: G {uj) ~ iG \u) = -iTrp 6{A 2 - e 2 ), 
with po = 1/(2A) (the real part of Go can be discarded 
to order A -1 ). This leads in the previous equations to 
integrals running from the lower bound e min = uj — A that 
can be replaced by e m i„ ~ — A for frequency smaller than 
A. This trick allows to obtain simple expressions for the 



self-energies: 

H f (u) = - 7 Po / de[G B (e)- J K ] + i 1Po [ deG' B {e) 

J-A J-A 

S s (w) = -po J "de G f (e) + ip J de G" f (e) (9) 



that lead by differentiation to the following relations: 

= ~7Po[G b (oj) - J K ] (10) 



d v r ,1 



-p Gf(aj) 



(11) 



Equation imply also the following boundary condition 
at the lower band edge (for T, B ): 



E B (-A)=ip f°deG'}(e) (12) 

J —A 



Following |l6j . we now define the inverse propagators 
Yf(uj) = —Gj(u) and Y b {uj) = G^j 1 (w). Combining the 
previous two relations with the Dyson equations (QJ-© 
for real frequencies, we easily arrive to: 



Id 1 - jpqJk d 

z-r — if — 





Yfdu f 



po doj B ^ Y B duj B ^ 



1 d 



Y f = ^T Y b 

J po ooj 



(14) 



The first equation provides a first constant of motion, 
that reads: 

log [ - Y f (u>)] +C= l ~ lpaJK Y B {uo) +1 \ogY B {uj) 

Po 

(15) 

where G is determined by the boundary condition 112fl . 
The crucial difference between the multichannel Kondo 
model we consider in this paper and previous work 0] 
on the NCA for the infinite U Anderson model is that, by 
inserting relation (|14fl into l|15|) . we obtain a first order 
differential equation: 



d_ 

duj 



Y n 



Poe C (r B r 7 e (1 ~ 7poJA " )yWpo 



(16) 



instead of a non-linear second order differential equation. 
Ultimately, this simplification is due to the disappearance 
of the high-energy scale associated with charge fluctua- 
tions on the impurity. We emphasize however that this 
result also applies to the usual NCA of the Anderson 
model, in the Kondo regime. 

To exhibit the exact solution of the multichannel 
Kondo model, we introduce two functions (x) defined 
by the integral (which is easily computed for integer 7): 



dt e ±4 V 



(17) 
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and thus obtain the following scaling form of the inverse 
bosonic propagator: 

Y b (uj) = Af«(oj/z K ) (18) 
where a — Sign(l — JPoJk) (19) 

The numerical prefactor in the previous equation is given 
by A = po/\l — jpoJk\, and a complex Kondo temper- 
ature zk = Pq 1 1 — 7A)^if I -1 " 7 e nas to be introduced. 
The fact that the propagators scale with respect to a com- 
plex energy scale stems from the presence of a spectral 
asymmetry for generic values of go jl2| • 

Equation l|18|) can be considered to be the complete 
solution of the large number of channels Kondo model 
for all energy scales below the high energy cut-off, and 
constitutes the central result of our paper. In particular, 
it is valid for frequency below and above the the Kondo 
temperature, defined as Tk = \%k\, as long as Tk is much 
smaller than the cut-off A. To analyze the properties of 
this scaling function is first particularly simple right at 
the Non Fermi Liquid fixed point. Indeed, equation Ijl6(l 
allows to locate the coupling J K = l/(po7) at which the 
various propagators exihibit scale invariance: 



Y B {oj) = U(l+7)e c w 



c, ,i V(i+7) 



at Jk = J 



K 



(20) 



for all frequencies smaller than the cut-off. We note also 
that, away from the fixed point J K , the behavior i|20|) 
still applies, albeit for uj < Tk- The power law ob- 
tained in equation (|2U[1 agrees with previous work [l2l[l6| . 
in which the zero-temperature low-frequency propagator 
G B {w) oc w 2Ab_1 , with 2A B = 7/(1 + 7), was obtained. 
For frequencies comparable to Tk and away from the 
fixed point, equation (|18fl leads to universal corrections 
to this pure power law behavior, a result that could 
not be obtained from the low frequency analysis per- 
formed in [l^. The various regimes observed are ex- 
pressed graphically in figure To understand first the 
universal crossover arising at Jk < Jk, we consider large 
frequency above Tk, and get Y B {u) oc log|w|, so that 
Gf(u>) oc 1/ui from equation i|14|) . This corresponds to a 
free moment regime, where the impurity spin is weakly 
bound to the screening cloud of the conduction electrons. 
The scaling function f~(x) thus describes the crossover 
from weak coupling to the fixed point J K [l|| ■ From the 
same token, f^(x) is associated to the flow from strong 
to intermediate coupling; we note however that in the 
case Jk > J* K one nas roughly Tk ~ A/r(7 +1), so 
that the condition Tk *C A (and hence universality of 
this regime) is met only for 7 > 1. This picture of the 
various crossover regimes that we have obtained is par- 
ticularly explicit, and illustrates the renormalization flow 
around the attractive fixed point. To make contact with 
the theory of quantum critical phenomena |l9j. we em- 
phasize that Jk is for the multichannel Kondo model 
an irrelevant perturbation around J K , so that no phase 
transition occurs as Jk is varied. Therefore the low en- 
ergy behavior at all Jk is controlled by the intermediate 



A 



Non universal 



Free Moment 



/ 



Unstable 
\ Over-Screened 
Moment 



NFL 



J 



J 



K 



K 



FIG. 1: Zero temperature (uj, Jk) crossover diagram: the 
dashed line represents the Kondo temperature Tk below 
which Non Fermi Liquid (NFL) behavior is observed; above 
this line two universal crossover regimes described in the text 
are characterized by the scaling functions 1181 . For \u)\ > A, 
non universal corrections start to appear. 



fixed point, in contrast to the case of a relevant tuning 
parameter driving a true quantum phase transition. 

Another way to depict the crossovers taking place 
around J K is to plot the scaling functions f^(x), fig- 



re n 





1 1 

/ 1 










/ 0.1 












0.01 0.1 


1 








1 



FIG. 2: Plot on both linear and logarithmic scales of the 
scaling functions (x) (lower curve) and /^T (x) (upper curve) 
for 7 = 1, as compared to their low energy asymptotic form 
[(1 + 7)z] 1/(1+7) (middle curve). 



It is also possible to compute the exact value for 
the Kondo temperature Tk- We will do this in the 
limit go — > 0> which corresponds actually to the usual 
Non Crossing Approximation. In this case, the pseudo- 
particle propagators get maximally asymmetric, and the 
spectral functions G b (uj) and G'j(u>) vanish for all uj < 0. 
This implies an explicit expression for the boundary con- 
dition l|12|) . T, B (— A) = 0, and a purely real inverse prop- 
agator Y b (lo) (for negative frequency). The initial con- 
dition Yg(— A) = 1/ Jk gives for Jk < Jk the following 
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relation: 



Tk 
A 



i/(poJk) — 1 



dt e* ti 



In the limit of small Jk , one finds: 

T K ~ A(p J^)T e T e - 1 /^ J ^ 



(21) 



(22) 



i.e. a vanishingly small Kondo temperature, as expected. 
When Jk approaches the fixed point value J K , Tk is 
maximum and Non Fermi Liquid properties apply up to 
the high energy cutoff. 

We now discuss temperature effects. In was shown 
how to obtain from the large- N equations I|2I5|) the exact 
finite-temperature low-frequency (i.e. to <C Tk) scaling 
functions using CFT arguments. It seems unfortunately 
impossible to apply this approach to extend our exact 
scaling function (fT%)) . valid for Tk ~ w <C A, to finite 
temperature. However we can succeed in generalizing 
the derivation of the linear differential equations (|10lll|) 
to the case of non zero temperature. For that purpose 
we consider the imaginary time saddle-point equations in 
the scaling regime: 



T, B {iv n ) 







-7 / dre tuJnT 

J (3 sin(7rr//3) 



n 



/3sin(7TT//3) 



[G b (t)-J k ] 
(23) 



This gives the idea of writing a finite difference equation 
for the self-energies, which provides indeed a remarkably 
simple result: 



T, f (iu n + in/ (3) - Y,j{iv n - inJJS) 

Sgfcjg + in /P) - T, B (iu) n - in/ (3) 
2n/f3 



= -i~fPo[G B {iv n ) - J K ] 
= -ip G f {iuj n ) (24) 



Although these relations do not lead to much analyti- 
cal understanding as compared to their zero temperature 
counterpart, they offer an economical way of tackling the 
finite temperature integral equations, which may be of 
interest for some generalizations of the Kondo model. 

We would like to conclude the paper on the possible 
applications of the formalism discussed in this work. Our 
main result, the computation of the full zero-temperature 
scaling functions of the many channel Kondo model, of- 
fers perspectives both in the physics of quantum impu- 
rities and strongly correlated systems. First, questions 
that we will address in future work concern the direct 
comparison of the exact scaling function l|18|l to the nu- 
merical solution of the large-N integral equations Q2I5JI . 
using both the standard numerical routines with fast 
Fourier fransforms and the system of finite difference 
equations l|24l) . In the same direction, it would also be in- 
teresting to perform a similar analysis and gain some an- 
alytical insight on the multichannel Kondo model in the 
bosonic representations |13| , a model that allows a richer 
phase diagram with a transition from over-screening to 
underscreening . 

As was pointed out in the introduction, models of sev- 
eral quantum impurities or lattice versions of the Kondo 
hamiltonian are of great experimental importance and 
notoriously difficult to tackle theoretically. We hope that 
the analytical step undertaken in the present work will 
allow further progress in this direction. 
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